Two-body problem in periodic potentials 
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We investigate the problem of two atoms interacting via a short range s-wave potential in the 
presence of a deep optical lattice of arbitrary dimension D. Using a tight binding approach, we derive 
analytical results for the properties of the bound state and the scattering amplitude. We show that 
the tunneling through the barriers induces a dimensional crossover from a confined regime at high 
energy to an anisotropic three dimensional regime at low energy. The critical value of the scattering 
length needed to form a two-body bound state shows a logaritmic dependence on the tunneling rate 
for D — 1 and a power law for D > 1. For the special case D — 1, we also compare our analytical 
predictions with exact numerics, finding remarkably good agreement. 
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I. INTRODUCTION 

Recent progress in controlling and manipulating the 
interatomic forces via Feshbach resonances |l| and the 
availability of tunable periodic potentials generated by 
laser beams are opening new fronts in the research on 
ultracold atomic gases. An important field of activities 
concerns the study of two-body physics in the presence of 
an external potential. This two-body problem is directly 
related to current experiments and moreover, it forms 
a basic step towards understanding many body states 
where strong correlations can play a significant role. 

From the experimental point of view, a major achieve- 
ment has been the formation of weakly bound diatomic 
molecules in experiments with two-component clouds of 
fermionic atoms 0- [J Q ■ Recently, also the binding 
energy of confinement induced molecules has been mea- 
sured by the group of Esslinger |7j working with samples 
of potassium 40 K in a deep two dimensional optical lat- 
tice. 

From the theoretical front, the two-body problem in 
the presence of harmonic trap ping , has been solved ana- 
lytically both for s-wa ve la, Of Eilini and p-wave short 
range interactions 0, Il3| . The case of a periodic ex- 
ternal potential is even more interesting as the center of 
mass and the relative motion no longer decouple |14lll5| . 
This gives rise to an interesting interplay between tun- 
neling and confinement effects. 

The properties of the bound state of two interacting 
atoms in the presence of a one dimensional optical lattice 
have been investigated numerically in Ref.[lJ|. For suf- 
ficiently deep lattices and binding energy Eb small com- 
pared to the interband gap, two different regimes were 
identified: a quasi two-dimensional regime for \Eb\ t, 
t being the interwell tunneling rate, where the molecular 
wavefunction is localized in one well and an anisotropic 
three dimensional regime for \Eb\ <C t where the wave- 
function spreads over many lattice sites |l6j . An analo- 
gous dimensional crossover occurs for the scattering am- 
plitude as the total energy of the two incident states be- 
comes large compared to the tunneling rate [l7j . 

In this work we apply the tight binding approach intro- 



duced in Ref.0 to study analytically the properties of 
the bound state and the scattering amplitude of two in- 
teracting atoms in the presence of a D dimensional tight 
optical lattice. We derive simple yet accurate predictions 
for many quantities that can be directly measured in cur- 
rent experiments with dilute ultracold gases. 

The paper is organised as follows. In Section II we 
present the general theory to solve the two-body problem 
in the presence of a periodic potential. In section III we 
apply the tight binding approach to the case D — 1 and 
we extensively compare our analytical results with the 
exact numerics of Ref . [l4| . In section IV we extend our 
analysis to the case of optical lattices of higher (D = 2, 3) 
dimension. 



II. GENERAL THEORY 

We consider two atoms in ordinary space interacting 
in the presence of a D dimensional periodic potential 
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where D = 1,2,3. Here s is the intensity of the laser 
beams generating the optical lattice, En = h 2 ir 2 /2md 2 
is the recoil energy, d is the lattice period and m the atom 
mass. There is no confinement in the remaining 3 — 1? 
directions. 

Modelling the interatomic interaction by a s-wave 
pseudopotential with coupling constant g = 4irh 2 a/m, 
with a the 3D scattering length, the Schrodinger equa- 
tion takes the form 
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where r = xi — X2 is the relative distance between the two 
atoms, R and ry are the components of the center of mass 
and the relative coordinates in the confined directions 
and V(R, ry) = V op t(x-i) + V op t(^2) is the total external 
field. 
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With respect to the case of harmonic trapping, the 
periodic field introduces a conceptually new difficulty re- 
lated to the fact that the center of mass and the rela- 
tive motion no longer separate. Nevertheless, the quasi- 
momentum Q, associated to the center of mass variable 
R, remains a conserved quantity even in the presence 
of interaction and can be used to classify the solutions 
of Eq.l©. This can be seen by noticing that Eq.JSJ) is 
invariant under the transformation 
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which shifts the center of mass of the two atoms by a lat- 
tice vector G leaving unchanged their relative distance. 



A. Bound state 

Bound states are solutions of Eq.J2J) whose energy E 
does not belong to the energy spectrum in the absence 
of interaction. In this case Eq.JSJ can be written in the 
integral form 

*(r,R) = J dR / G B (r J R;0 J R , )fl^ 7 (r / *(r , ,R , )) r , = o, 

(5) 

where Ge is the Green function associated to Eq. (J2J with 
g = 0. The behaviour of the Green function at short dis- 
tance is dominated by the kinetic term in Eq. J2J ■ Taking 
into account that for small r 



el Q(R-R') e *kr d 3 k d DQ 

D 



h 2 k 2 /m + h 2 Q 2 /4m (2tt) 3 (2tt) 
r— r') d D Q 



A-rth 2 ! 



iQ( 



(2tt) b ' 



(6) 



the Green function admits the following expansion 



equation. Second, the conservation of quasi-momentum 
in periodic potentials is ensured by the symmetry prop- 
erty K E (R,R!) = K E (R + G,R + G) of the kernel, 
holding for any lattice vector G. 

In the following we are mainly interested in the weakly 
bound state whose energy E is below the minimum E re t 
of the non interacting energy spectrum 

Hi- 

For a given 

quasi-momentum Q, the lowest energy solution of Eq.J2J) 
for g — corresponds to the case where both atoms oc- 
cupy the state of quasi-momentum Q/2 in the lowest 
Bloch band, that is E ref (Q) = 2e x (Q/2). We define the 
binding energy of the molecule as Ef, = E — E re f < 0. 

A direct consequence of the non separability of the 
center of mass and the relative motion in the lattice is 
that the binding energy becomes a function of the quasi- 
momentum of the molecule, or, the other way around, the 
tunneling properties of the molecule (strongly) depend on 
the value of the scattering length. A second important 
difference from the case of harmonic trapping is that, in 
the lattice, Eq.@ predicts a finite critical value of the 
scattering length a = a cr < such that no bound state 
exists for a cr < a < 0. Notice that in free space (s = 0), a 
bound state exists only for positive value of the scattering 
length, i.e. a cr — — oo. 



B. Scattering amplitude 

Let us consider a non interacting two-particel state 
* (r,R) = 0n 1 ,q 1 (ri)0n 2 ,q 2 (r 2 )e lkiri , where n , q (x) = 



-/ (R) = / dR! K E (R,R')f(R'). 
9 J 



(8) 



The solution of Eq.© yields the energy E as a function 
of the 3D scattering length a. 

Two comments are in order here. First, the integral 
equation (JSJ is completely general and is valid for any 
external potential. In the special case of harmonic trap- 
ping, /(R) is also an eigenstate of the center of mass 
Hamiltonian and Eq.Q reduces to a simpler algebraic 



T\f=i u nj,qj(' x )- Here u n ^ q {z) are the solutions of the one 
dimensional Hamiltonian H = —(ti 2 /2m)d 2 /dz 2 + V op t(z) 
with energy e n (q), with n the band index and q the quasi- 
momentum. A quantity of great physical interest is the 
scattering amplitude associated to the incident state. 
In the presence of the lattice, we define it as 

/ sc [*o]=a y**(0,R)a r (r*(r,R)) r=0 dR, (9) 

where a is the 3D scattering length and \t(r, R) is the 
solution of Eq. , which can be written as 

¥(r,R) = # (r,R) (10) 

+ J dR , G E+l o(r, R; 0, R')3^7 (r'*(r', R')W ■ 

Here G E +io is the retarded Green function evaluated 
at the energy E — e ni (qi) + e„ 2 (q2) + h 2 k\/m of the 
incident state and e n (q) = J2j=i e «j (ij)- Introducing 
/(R) = d r (r\&(r, R)) r=0 and making use of the expan- 
sion J7J) in Ea. i|10[l . we find the inhomogeneous integral 
equation 

/(R) = *o(R,0) + 5 j dM K E+i0 (K,Rf)f(R'). 

(11) 

which has to be solved to find the scattering amplitude 
@. The comments below Eq.© applies here as well. 



G E (r, R; 0, R) = —^-S D (R - R')+K E (R, R')+0(r), 

where K E (R, R') is a regular kernel which depends on 
the energy and the external potential. 

When inserted into Eq.©, Eq.(J7J| yields the Bethe- 
Peierls boundary condition $(r, R) ~ (a/r — 1)/(R) for 
r — * 0, where /(R) is a function of the center of mass 
position satisfying the integral equation 
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C. Numerical and tight binding solution 

Equations J3J) and (|llfl can be solved numerically fol- 
lowing the method developed in Ref.^iJ for the case of a 
ID optical lattice. In order to extract the regular kernel 
Ke from the Green function, the latter is expanded in 
the basis of the non interacting states 

G £ (r,R;0,R') = (12) 
d 3 ~ D k ± p d°^_ * k± . r± 



E 



ni,n 2 
0ni,qi(R 



(27T) 



q2 (R-^)< 1)qi (R')^ 2!q2 (R') 



E - e ni (qi) - e n2 (q 2 ) - h 2 k\/m 

Since the singular term in the rhs of Eq.JJJ) docs 
not depend on the external potential, we add and 
substract from Ge, the Green function l|12(l eval- 
uated for s — 0. This permits us to write 
the kernel as K E (R,K') = lim r ^o[G_e(r, R; 0, R') - 
G|r (r,R;0,R')] + K S E =°(R, R), where the latter is 
given by 



Kff°(R,R) 



dP 



4irh 2 J (2tt 



e lP(R - R,) y/\E\ + h 2 P 2 /4n 

The term in square brackets is evaluated numerically 
from Ea. l)12|l . Once the kernel K E is known, Eqs (jHJ) 
and I f! IP can be solved using standard diagonalization 
routines. This numerical method has been applied in 
Ref. f° r the case D = 1. It should be noticed that for 
D > 1 the numerical effort needed to solve Eqs © and 
jll) , using the above method, grows considerably. 

In this paper we restrict ourselves to the case of suf- 
ficiently deep lattices and energy E small compared to 
the interband gap e g . In this limit the energy depen- 
dence of the Green function in Ea. (|12|) comes entirely 
from the lowest Bloch band and the two-body problem 
can be solved perturbatively using a tight binding ap- 
proach pioneered in Ref. |l7j . 



III. ONE DIMENSIONAL LATTICE 

Below, we discuss in details the case of a one dimen- 
sional optical lattice (D = 1) by comparing our analytical 
predictions with the exact numerics of Ref. |14| . 

A. Bound State 

We start our analysis by noticing that for D = 1 the 
integration over kj_ in Eq. (|12fl is ultraviolet divergent for 
rj_ — > 0. Taking into account that the singular term in 
the rhs of Eq. J7J) does not depend on energy, we can write 
identically: 

K E (Z,Z') = [G E (r,Z;0,Z') - G E=Eref (r,Z;0,Z% =0 
+ K E=Eref (Z,Z'). (14) 



We consider first the term between square brackets which 
contains the energy dependent part of the kernel. The 
integration over kj_ in Eci. (|12|) converges even for rj_ = 0. 
Moreover, for energy E small compared to the inter- 
band gap, only the lowest band n\ — rii — 1 con- 
tribute significantly. For sufficiently large values of the 
laser intensity, the states of the lowest Bloch band can 
be written in terms of Wannier functions as <j>i qz (z) ~ 
e llq ^ d w{z — id), where 



w(z) = 



1 



r 1 /4 (T i/2 



is a variational gaussian ansatz. By minimizing the en- 
ergy J dzw{z)[—(h 2 /2m)d 2 /dz 2 + V opt (z)]w(z) with re- 
spect to a, one finds d/a ~ its 1 / 4 exp(— 1/Ay/s). 

Substituting the tight binding expression for the Bloch 
states in Ea. H12|) . we immediately see that the term in 
square brackets is diagonalized by the ansatz 



/ Q (Z)~]>> 2 (Z-^)e 



iQZ 



(16) 



where Q is the quasi-momentum of the molecule. Moti- 
vated by this fact, we have verified numerically that the 
ansatz <(TT)|) is a solution of Eq.(|HJ with E — E re f, that 
is 

J dZ'K E = Eref (Z, Z')f Q {Z')dZ' = f Q (Z)m/4irh 2 a cr , 

where a cr = a cr (Q) is the critical value of the scatter- 
ing length needed to form a two-body bound state with 
quasi-momentum Q. The ratio d/a cr as a function of the 
laser intensity has been evaluated numerically in Ref. 0] 
and in Ref. |17| | an analytical estimate was derived for the 
special case Q = 0. 

By inserting the ansatz l|16|l into Eq. JSJ and integrating 
over Z' , we obtain 



— J V2na = -a 1D (E b , Q), 



a a c 



where 



a\D{Eb, Q) = 



^ln 
2tt 



-E b -2e 1 (Q/2)+e 1 (q)+e 1 (Q- 



(18) 



(19) 



ei(g)+ei(Q-?) 



and we have used the definition of binding energy E b = 
E — 2ei(Q/2). Equations i|18ll9fl give Eb as a function of 
the quasi-momentum Q, the scattering length a and the 
lattice parameters. 

Let us first discuss the case of zero quasi-momentum. 
By substituting the tight binding dispersion t\{q) = 
2t(\ — cos(qd)) in Ea. (|19|l and integrating over q, we ob- 
tain 



a 1D (E b ,0) 




(20) 
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FIG. 1: Binding energy as a function of the inverse scattering 
length for different values of the laser intensity: the analytical 
prediction (solid line) is compared with the exact numerical 
data (squares) . Also shown are the position of the bandwidth 
(dashed line) and the asymptotic behaviours for the binding 
energy (dotted lines) in the two regimes [see Eqs 12111 and 



Equation l|20|) describes the dimensional crossover in the 
binding energy as the ratio \E b \/At is varied. The three 
dimensional regime \E b \ <C 4i is charachterized by a 
molecular wavefunction extended over many lattice sites. 
Expanding Eq. l|20ll around E b = and inserted it into 
Eq.l^HJ, we nna 



1/a - l/a cr = C^/\E b \m*/h, (21) 

where m* = h 2 /2td 2 is the atomic effective mass eval- 
uated at the bottom of the band and C — d/^/2n<T. 
The above result shows that the optical lattice gives rise 
to an effective shift of the resonance from 1/a = to 
1/a = l/a cr < 0. In the quasi 2D regime 4t <C \Eb\ -C e g , 
the two interacting atoms are localized at the bottom of 
the same optical well, where, to a first approximation, the 
potential Q is harmonic with frequency uiq = h/ma 2 . In 
this case, the center of mass and the relative motion de- 
couple and the problem can been solved analytically. In 
the relevant limit \E b \ *C ujq ~ e ff , the binding energy is 
given by 

E%° = -Aftw exp(-v / 2^CT/|a|) (22) 

with A = 0.915/tt 10]. In Fig Q] we plot the binding 
energy versus inverse scattering length derived from Eqs 
(|18fl and (|20|l for different values of the laser intensity 
(the value of the parameter d/a cr has been taken from 
Ref.^H). We find that the analytical predictions (solid 
line) are in perfect agreement with the numerical data 
(squares) provided the binding energy is small compared 
to the interband gap. This condition is no longer satisfied 
approaching the resonance from the negative side, so our 
method is accurate only for scattering length a < 0. 



FIG. 2: Critical value of the inverse scattering length as a 
function of the laser intensity for quasi-momentum Q — 
and Q = qb- The analytical prediction (solid line) is shown 
with the exact numerical data (circles). 



Interestingly enough, the fact that E b ~ E> b w for 
Ei, 3> it can be used to derive an analytical expression 
for the critical scattering length . By substituting the 
expansion am(Eb, 0) ~ -ha(\E b \/2t) in Eq.JTSJ with E b 
given by Eq.^, we find [jj 

±(Q = 0)=*ln(^). (23) 
a cr v27rcr \ tt J 

This is plotted in FigEl as a function of the laser in- 
tensity (solid line) together with the numerical results of 
Ref.[14| (circles). We see that the agreement is excellent 
already at relatively low values of the laser intensity. The 
approximate linear behaviour as a function of s shown in 
FigEl is due to the logarithmic dependence of the ratio 
(|23|l on the tunneling rate t and is a specific feature of 
the case D = 1. It is worth pointing out that a good 
estimate of be obtained by the following physical 

argument. In the presence of a deep lattice, the binding 
energy is basically given by Eq. l|22|l provided \E b \ ^> At. 
The tunneling effects that cause the molecule to break 
become important when the binding energy (in modu- 
lus) is of the order of the bandwidth 4t or smaller. This 
suggest that the critical scattering length a = a cr can 
be estimated from Eq.O by setting —E^° ~ At. This 
reproduces the result (|23|l in the limit t <C Hujq. 

The effective mass M* of the molecule is defined as 
h 2 /M* = d 2 E(Q)/0Q 2 evaluated at Q = 0. From Eqs 
(|18I19|) after some algebra we obtain 

2m * i x 1^1 l(E b \ 2 \E b \ 

-NF = l + ^r-i\-t) (24) 

where E b = E b (0) and m* is the effective mass for single 
atoms. We see that in the anisotropic three dimensional 
regime, the mass ratio 2m* /M* = 1 for E b = and de- 
creases as —E b increases. In the quasi two dimensional 
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FIG. 3: The inverse effective mass of the molecule as a func- 
tion of the inverse scattering length for s — 5: analytical 
prediction (solid line) and exact numerical data (squares). 




d/a 



FIG. 4: Molecular bandwidth as a function of the inverse 
scattering length for different values of the laser intensity: 
analytical prediction (solid line) and exact numerical data 
(squares). The atomic bandwidth is w at = 0.264Er for s — 5 
and Wat = Q.077E R for s = 10. 



regime, -E b > At, Eq.J23 yields 2m* /M* = 2t/\E b \ < 
1, showing that the correlation between the two atoms 
results in an increased inertia of the molecule. The mass 
ratio i|24|) is plotted in Fig|3]against the inverse scattering 
length for s = 5 togheter with the exact numerical result 
(dots). We emphasize that the dependence of the effec- 
tive mass of the molecule on the value of the scattering 
length is a clear consequence of the non separability of 
the center of mass and the relative motion in the lattice. 

Let us now discuss the solution of Eq. J5J) corresponding 
to finite value of the quasi-momentum Q. The binding 
energy Eb(Q) can be obtained from Eqs l|18|) and i|19|) 
provided the critical scattering length a cr {Q) is known. 
This quantity can be easily calculated by noticing that 
in the quasi two dimensional limit \Eb(Q)\ ^> 4t, to low- 
est order in tunneling, the binding energy does not de- 
pend on the center of mass motion and therefore Eb(Q) = 
Eb(0) = Elf . We have found that for Q ~ q B the use 
of the tight-binding dispersion e\(q) = 2t(l — cos(qd)) in 
Eq. (|19JI gives inaccurate results. This dispersion relation 
predicts indeed a continuous degeneracy for the lowest 
energy state with quasi-momentum Q = q B of the two 
atoms in the absence of interaction. As a result, the in- 
tegration over quasi-momentum in Eq. i|19fl is divergent 
at this point. This artefact can be easily eliminated by 
using the numerically exact dispersion or by including 
higher harmonics. With this important remark, the cal- 
culated critical value of the scattering length for Q = qs 
agrees fairly well with numerics, see Fig|3 

Another quantity of physical interest is the molecular 
bandwidth w = E(Q = q B ) — E(0). Within our notation, 
the latter can be written as w — E re f(q B ) + Eb(q B ) — 
Eb(0). This is plotted in Fig0]as a function of the inverse 
scattering length. Similarly to the effective mass, the 
bandwith decreases as d/a increases. However, for s = 5 
(left panel) there is a discrepancy between the analytical 
prediction and the numerical result which was not found 
in Fig|S] This is due to the fact that the interband gap 
is smaller at the edge of the Brillouin zone and there 



is a residual contribution to Eb(qB) coming from higher 
bands. This effect disappears when s is increased so that 
for s — 10 good agreement with the numerical results 
is obtained (right pannel). Notice that, the molecular 
bandwidth is always smaller than the atomic bandwidth 
w at . This can be seen from the fact that Eb(q B )—Eb(0) < 
and E ref (q B ) = 2e 1 (q B /2) < ei(q B ) = w at . 



B. Scattering amplitude 

We assume that the two scattering atoms belong to 
the lowest Bloch band, so the incident state in Eq. I|llfl is 
given by ^ (r,Z) = <j> lqi (zi)0i 92 (z 2 )e* k±r± . The energy 
is E = h 2 k'j_/m+ei(qi)+ei(q2) and the quasi-momentum 
Q is given by 



Q = 9i + 92 if qi 
= qi + q 2 ± 2q B if 



- <?2 < QB 

\qi+q2\>qB- (25) 



Notice that we only require E <C e ff , but there is no 
restriction on the value of the scattering length a. We 
insert the ansatz f Q (Z) = AY /t w 2 (Z - £d)e lQZ in 
Eci. (|ll|) where A is a numerical coefficient which is re- 
lated to the scattering amplitude © by f sc = aCA, 
where C = dj \f2~Ko. Repeating the same steps as in 
the previous section, after integration over Z' ', we find 



fsc - 

where 
Pid(E,Q) 



aC 



1 - a/a cr + {a/^/2-Ka)f3 1D (E, Q) 



(26) 



T dq 
r * dq 



hi 



-E + ei(q) + ei(Q - q) 



ei(g) + d(Q-g) 



m j_ ^0(£-ei(<z)-ei(Q-e?)).(27) 
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This function is related to the analytical continuation of 
ot\r>i defined in Ea. (|19|l . Here 8(x) is the unit-step func- 
tion and we have used the formula (x — E— i0) _1 = P(a; — 
E)^ 1 + m8(x — E), where the symbol P stands for Princi- 
pal value. In the special case Q — 0, the integration over 
quasi- momentum in Eq. (|27(l can be done analytically us- 
ing the dispersion relation ei(q z ) — 2t(l — cos(q z d)) and 
one finds 17] 

f3 1D (E, 0) = i arccos(l - E/tt) if E < 8t, (28) 
- ]n[E(l + y/l - 8t/E) 2 /8t] +iTr ifE> 8t. 

Equations (|27|l and l|28|) show that the scattering am- 
plitude undergoes a dimensional crossover as the ra- 
tio E/8t is changed. In the anisotropic three dimen- 
sional regime E <C 8t, one has f sc = aC/(l — a/a cr + 
iaCy/Em* /h). Taking into account Eq. (|23|) . in the 
quasi two dimensional regime E 3> 8t we obtain f sc — 
d(a/V^<y)/ (l + (a/v / 2Tc^)(ln[ASa;o/£ , ] + wr)) and A = 
0.915/-7T in agreement with Ref.[To|. 

It is interesting to notice that the scattering amplitude 
for Q = is related to the effective coupling constant g e f / 
for Fermi atoms undergoing Cooper pairing in the case of 
weak attractive interaction (a < 0). This result was used 
in Ref. 17] to calculate the superfluid transition temper- 
ature of a two-component Fermi gas in the presence of a 
one dimensional optical lattice. 

IV. TWO AND THREE DIMENSIONAL 
LATTICES 

In this section we extend our analysis to higher di- 
mensional lattices. Differently from the case D = 1, the 



integration over k± in Eg. I| 12(1 for D = 2,3 is convergent 
for rj_ — 0, which considerably simplifies our analysis. 
Moreover, in the limit E <C fuvo, the energy dependence 
of the kernel comes entirely from the contribution of the 
lowest band (ni — ri2 — 1). The inclusion of higher 
bands gives small correction to the binding energy but is 
crucial to reproduce the correct behaviour of the scatter- 
ing amplitude close to Feshbach resonance, as discussed 
below. 



A. Bound state 

For a D dimensional optical lattice, the tight binding 
ansatz (I16fl for the function /(R) in Eq.JSJ generalizes 
to 



f(R)=UP =1 f Ql (R i ), (29) 

where Qi and Ri are the «-th components of the quasi- 
momentum Q and the center of mass position R, re- 
spectively. By inserting the ansatz (|29|l in Eq.JSJ and 
integrating over R', we find 

l/g = A D (E b ,Q) + B D (Q), (30) 

where Ad gives the contribution of the lowest Bloch band 



A 2 (£b,Q) 
A 3 (# b ,Q) 



-C 3 



'IB 



d 2 q 



1 



-9b ( 27r ) 2 V- E b - 2ei(Q/2) + ei(q) + ci(Q - q) ' 

qB jPq 1 

- qB (2tt)3 -E b - 2 Cl (Q/2) + ei (q) + ei(Q - q) 

I 



(31) 
(32) 



and C has been defined previously. The quantity Bd(Q) 
in Ea. H30(l takes into account the small correction com- 
ing from the contribution of all higher bands and can be 
evaluated cither numerically or analytically, as explained 
below. Eqs (|30I32|) permit us to calculate the binding 
energy as a function of the scattering length for a given 
quasi- momentum Q. 

For simplicity we restrict our discussion to the case 
Q = 0. In analogy with the case D = 1, there are 
two regimes: an anisotropic three-dimensional regime for 
£jCi and a confined regime of reduced dimensionality 
for t <C \Eb\ <C huiQ r~j e The constant Bo can be fixed 



by matching the binding energy in the limit \E b \ t with 
the corresponding result valid for the harmonic oscillator 
potential 8, 9] in the limit \E b \ <§C ujq- 



\E h 



E h = —— \ — 



m a i (i 

h 2 [2 



V2 



in 



7T O" 



(1- 



m-, 



(D=2) 
(D=3), 



(33) 
(34) 



where r] 2 = 1.0326 and 773 = 0.2448. This gives B D = 
•qoTn/Anh^a. Taking into account that the range of ap- 
plicability of our method is \Eb\ "C hwo or, equivalently, 
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\a\ <C cr, we see that the contribution from higher bands 
to the binding energy is typically small. 

The critical value of the scattering length a — a cr can 
be obtained by setting E b = in Eqs (|3lT) and 
and integrating over the quasi-momcnta using the tight- 



binding dispersion e n (q) 
gives 



Ef=i2t(l 



cos(qjd)). This 



(D = 2) (35) 




where I2 = 0.454 and I3 = 0.101. Similar results for 
D = 3 have been obtained by Fedichev et al.|lj| using a 
different method. Similarly to the ID case, the functional 
dependence of the critical scattering length on the ratio 
hu /t simply follows from Eqs [D=2] and (JSU [D=3] 
by substituting the binding energy with the tunneling 
rate (\E b \ — > t). We see from Eqs dSSJ and ^ that the 
ratio d/a cr has a power law dependence on the tunneling 
rate t and therefore it increases exponentially fast as a 
function of the laser intensity, in contrast with the case 
D = 1 shown in Fig|5] 

From Eqs l|31Jl and (|32|l . we recover the behaviour of 
the binding energy close to the critical point 



1 

a 



1 



*D/2 
7j m (D-l)/2' 



which characterizes the three dimensional regime. 

The effective mass M* for the molecule (at Q = 0) can 
be calculated from Eqs (|3T|l and ijlfffy. For zero binding 
energy, the mass ratio 2m* /M* = 1 and decreases rapidly 
by increasing \E b \ as an effect of the correlated motion of 
the two atoms. In the confined regime \E b \ 3> t, we find 
2m*/M* = (D + l)t/\E b \. 



B. Scattering amplitude 



We assume again that the two scattering atoms belong 
to the lowest Bloch band, so the incident state in Eq. (jl 1|> 
is given by * (r,R.) = <?W (xi)<?W {x 2 )e lk±T± and its 
energy is E = h 2 k 2 _/m+e 1 (qi) + e 1 (q2)- By inserting the 
ansazt /q(R) = AH^L x/q^-RO, in Ea. ((TTJl and taking 
into account that f sc — aC D A, we find 



JSC — 



aC L 



l-g(p D (E,Q) + B D y 



(38) 



(37) where Bp has been defined above and 



02(E,Q) = 



-C 2 ^ f qB d 2 q 



2h 



-c* 



'ID 



IB 

d 3 q 



9 (ei(q) + ei(Q - q) - E) ? 9 (E - ei(q) - d(Q - <Q) 
1 



V^-ei(q)-ei(Q-q) 



, s (27r)3-S- l 0+(Q/2) + e 1 (q) 

I 



£i(Q 



(39) 
(40) 



Equations l|38|) and H40|l permit us to calculate the low en- 
ergy scattering amplitude f sc = aC 2 A in the presence of 
the lattice. Notice that f sc has a pole at the critical value 
of the scattering length a — a cr (Q) for E = £ , re /(Q), 
corresponding to the bound state with vanishing binding 
energy. 

For simplicity we confine ourselves to the case Q = 
where E re f = 0. In the three dimensional regime, from 
Eqs iffify and we find the typical result 



JSC 



aC 



D 



(41) 



1 - a/a cr + iaC D m* D / 2 VE/m( D - i y 2 h 

valid in the limit E <C t. This result was also found 
in Ref. ^(| for the case of a 3D lattice. In the case 
D = 2 the system undergoes for E > 8t a crossover to a 
quasi-one dimensional regime. The dispersion of the low- 
est band can be neglected in Eq. l|39[l and the scattering 



amplitude © reduces to the asymptotic value 

aC 2 VE 



E + i 



(42) 



\frna' 2 



The inclusion of higher bands (772 7^ 0) becomes crucial 
when the scattering length is large and only when they 
are taken into account, one recovers the confinement in- 
duced resonance for a — ojr\2- Notice that our defini- 
tion of the scattering amplitude [see Eq. ©] is related to 
the amplitude f(E) used in Ref.@ by a scaling factor 
f(E) - -i{V2^/d 2 VE^)f sc {E). 

Finally, we have verified that for D > 2, Eqns (|30|l and 
(|38|l also follow from a Hubbard-like model with hopping 
along the the optical lattice, free motion in the perpendic- 
ular direction and the contact interaction strength equal 
to U = £//[(! - 7] D a/a)(V^cr) D ] HJ. 
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V. CONCLUSIONS 

In this paper we have generalized the two-body scat- 
tering theory to include the effect of periodic fields, where 
the center of mass and the relative motion do not sep- 
arate. The interplay between confinement and tunnel- 
ing effects gives rise to an interesting crossover between 
a low dimensional regime where the interwell hopping 
can be treated perturbatively and an anisotropic three 
dimensional regime where the wave function for the rel- 
ative motion spreads over many lattice sites. By using 
a tight-binding approach, we have investigated the low 
energy properties of the bound state and the scattering 



amplitude across the two regimes. The precision of our 
analytical predictions have been tested for the case of a 
one dimensional lattice where exact numerics is available, 
finding very good agreement. Our results are relevant for 
current experiments on ultracold gases in optical lattices. 
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